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THÉORIE DES FONCTIONS NUMÉRIQUES SIMPLEMENT 

PÉRIODIQUES. 

Par Edouard Lucas, Professeur au Lycée Cliarlemagne, Paris. 

Ce mémoire a pour objet l'étude des fonctions symétriques des racines 
d'une équation du second degré, et son application à la théorie des nombres 
premiers. Nous indiquons dès le commencement, l'analogie complète de ces 
fonctions symétriques avec les fonctions circulaires et hyperboliques ; nous 
montrons ensuite la liaison qui existe entre ces fonctions symétriques et les 
théories des déterminants, des combinaisons, des fractions continues, de la 
divisibilité, des diviseurs quadratiques, des radicaux continus, de la division 
de la circonférence, de l'analyse indéterminée du second degré, des résidus 
quadratiques, de la décomposition des grands nombres en facteurs premiers, 
etc. Cette méthode est le point de départ d'une étude plus complète, des pro- 
priétés des fonctions symétriques des racines d'une équation algébrique, de 
degré quelconque, à coefficients commensurables, dans leurs rapports avec les 
théories des fonctions elliptiques et abéliennes, des résidus potentiels, et de 
l'analyse indéterminée des degrés supérieurs. 

Section I. 

Définition des fonctions numériques simplement périodiques. 

Désignons par a et b les deux racines de l'équation 

(1) & = Px — Q , 

dont les coefficients P et Q sont des nombres entiers, positifs ou négatifs, et 
premiers entre eux. On a 

a + b = P , ah — Q; 

et, en désignant par S la différence a — b des racines, et par A le carré de 
cette différence, on a encore 

_P + è x_-P-3 



a = t_J^ , b = ^—^ , 8 = VAz=VP 2 — 4Q. 
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Cela posé, nous considérerons les deux fonctions numériques JJ et F définies 
par les égalités 

(2) U n = a ° T , V n = a»+b». 

a — o 

Ces fonctions U n et V n donnent naissance, pour toutes les valeurs entières 
et positives de n, à trois séries d'espèces différentes, selon la nature des racines 
« et 6 de l'équation (1). Cette équation peut avoir : 

1°. Les racines réelles et entières ; 

2°. Les racines réelles et incommensurables ; 

3°. Les racines imaginaires. 

Les fonctions numériques de première espèce correspondent à toutes les 
valeurs entières de a et de b, et peuvent être calculées directement, pour 
toutes les valeurs entières et positives de n, par l'emploi des formules (2). 
Si l'on suppose plus particulièrement a = 2 et b zz 1, on trouve, en formant 
les valeurs de U„ et de V„, les séries récurrentes 



ni 



n: 0,1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, , 

U n : 0, 1, 3, 7, 15, 31, 63, 127, 255, 511, 1023, 2047, , 

V n : 2, 3, 5, 9, 17, 33, 65, 129, 257, 513, 1025, 2049, 

étudiées pour la première fois par l'illustre Fermât. Nous observerons, 
dès maintenant, que la série des V n est contenue, pour les trois cas que nous 
considérons, dans la série des U n , puisque les formules (2) nous donnent la 
relation générale 

(3) ' U 2n =U n V n . 



Les fonctions numériques de seconde espèce correspondent à toutes les va- 
leurs incommensurables de a et de b dont la somme et le produit sont 
commensurables. On peut les calculer en fonction de la somme P et du 
discriminant A de l'équation proposée, au moyen des formules suivantes. Le 
développement du binôme nous donne 

1 ï . 2 1.2.3 " " ' ' 

2»M P» n pn — ljï _t_ n { n -V pn — 2(j2 n V 1 £j V 1 " ) P»-3M I ( ^\» . 

- ~ x 0_r 1.2 1.2.3 •••TV o;, 

et, par soustraction et par addition, 

47 
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(4) 



9»-itt _ n p.-i.i, n(n — l)(n — 2) pm _ 3A 

w ( W _l)( w _2)(w — 3) (m — 4) pK _ s 
+ 1.2.3.4~5 A + . . . . , 

g.-* Fm = pn + nJn-Vl pn _ 2A + n(n-l)(n--2)(n-3) p^ + _ 
1.2 1 . 2 . o . 4 

On obtient ainsi, pour les premiers termes, 

C = 0, tfizrl, Z7 2 = P, U a = P*—Q, U i = P* — 2PQ, 

F =2, F 1= P, F 2 = P 2 — 2Q, V 3 = P S — 3PQ, F 4 = P 4 — 4P 2 Q + 2Q*. 

Les fonctions numériques de seconde espèce les plus simples correspondent 

aux hypothèses 

P=l, Q — — 1, A = 5, 

ou à l'équation a? = oc + 1 ; 

on a, dans ce cas, 

.37t 1 + V 5 , o- 71 1 — VIT 

a = 2 sm — - = — -pr , o = — 2 sin — = , 

10 2 10 2 

et, par suite, en désignant par u n et v n les fonctions qui en résultent, 

(1 + V5)»— (1 - V^) M .,_(! + VJ3) re + (1 - V 5) " 

On forme ainsi, pour les premières valeurs de n entières et positives, les 
séries 

n: 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, , 

u n : 0, 1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, , 

v n : 2, 1, 3, 4, 7, 11, 18, 29, 47, 76, 123, 199, 

La série des u n a été considérée pour la première fois par Léonard Fibonacci, 
de Pise.* Elle a été étudiée par Albert Gtirard,| qui a observé que les 
trois nombres u m u M u n + 1 forment un triangle isoscèle dont l'angle au sommet 
est à fort peu près égal à l'angle du pentagone régulier. Robert Simson J a 



* Il liber Abbaci di Leonardo Pisano, pubblicato seconda la lezione del Codice Magliabechiano, da B. Boncom- 
pagni. Koma, 1867. Pag. 283 et 284. 

f L'arithmétique de Simon Stevin, de Bruges, revue, corrigée et augmentée de plusieurs traictez et annotations 
par Albert Girard, etc. Leide, 1633. Pag. 169 et 170. 

J Philosophical Transactions qf the Royal Society qf London, Vol. xlviii, Part I, for the year 1753. An 
explication of an obscure passage in Albert Girard's Oommentary upon Simon Stevin's Works. Pag. 368 
et suiv. 
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fait remarquer en 1753, que cette série est donnée par le calcul des quotients 
et des fractions convergentes des expressions irrationnelles 

V~5 + 1 . VI- 1 

! — et . 

2 2 

En 1843, J. Binet * donne, au moyen de cette série, l'expression du dénom- 
brement des combinaisons discontigues. En 1844, LAMÉf indique l'appli- 
cation que l'on peut faire de cette série à la détermination d'une limite 
supérieure du nombre des opérations à faire dans la recherche du plus grand 
commun diviseur de deux nombres entiers. 

Nous prendrons aussi quelquefois pour exemple la série U n de seconde 
espèce, donnée par les hypothèses 

A = 2 2 .2, 



P=2, 




Q = — 1, 


ou par l'équation 




x 2 = 2x + 1. 


On a alors les séries 






n: 0,1, 2, 3, 4, 


5, 


6,. 7, 


U n : 0, 1, 2, 5, 12, 


29, 


70, 169, 


V n : 2, 2, 6, 14, 34, 


82, 


198, 478, 



8, 9, 10, 11, . . . 

408, 985, 2378, 5741, . . . 
478, 1154, 2786, 6726, 16238, . . . 

que nous désignerons sous le nom de Séries de Pell, en l'honneur du 
mathématicien de ce nom qui résolut, le premier, un célèbre problème 
d'analyse indéterminée proposé par Fermât, et concernant la résolution 
en nombres entiers, de l'équation indéterminée 

x 2 — Af = ± 1. 

Les fonctions numériques de troisième espèce correspondent à toutes les 

valeurs imaginaires de a et de & dont la somme et le produit sont réels et 

commensurables. Les plus simples proviennent des hypothèses 

P-l, Q = l, A = —3; 

on a, dans ce cas, 

1 + v^="3 . 1 — V^3 
« = , b = _ , 

par conséquent a et b sont les racines cubiques imaginaires de l'unité négative ; 
de plus, 

* Comptes rendus de l'académie des sciences de Paris, tome, xvii, pag 562 ; tome xix, pag. 939. 
f Comptes rendus, etc., tome xix, pag. 867. 



188 Lucas, Théorie des Fonctions Numériques Simplement Périodiques. 

Ainsi les valeurs de U n reviennent périodiquement clans l'ordre 

0, 1, 1, 0, -1,-1, 

et donnent lieu à un grand nombre de formules simples déduites des pro- 
priétés générales des fonctions U n et V n , et concernant la trisection de la 
circonférence. 

Quelquefois aussi nous considérerons les séries analogues déduites de 

l'équation 

x 2 — 2x — 2, 

dans laquelle 

a = 1 + V^T, 5 = 1 — V — 1, A = — 2 2 , 

et les séries déduites de l'équation 

x 2 = 2x — 3, 
dans laquelle 

a = 1 + V"^2, b = l — V^2, A = — 2 X 2 2 ; 

nous désignerons les séries obtenues dans cette dernière hypothèse, sous le 
nom de séries conjuguées de Pell. 



Section II. 
Des relations des fonctions U n et V n avec les fonctions circulaires et hyperboliques. 
Si l'on fait 



n T a 

z = — - Los-, nep. — 

2 ° l b 



dans les formules 



cos 



(z V — 1) 



_e° + e~ 



sin (z V — 1) .= 



2 
— e~ 



on obtient 



S1D 



Log. £) = 



2 V — 1 
1 



n *~ n 

Lb* a\ 



cos (nV-J) 

V 2 

/ n V — 1 T a \ 1 

(—2— L ° g - T) = 2V=Ï 






n n 
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on a donc, entre les fonctions U n et V m et les fonctions circulaires, les deux 
relations 



(5) 



V n = 2Q* cos ( — Log. — ) , 



TT 2(f . (nsj — 1 T a\ 



Il résulte immédiatement de ce rapprochement que chacune des formules de 
la trigonométrie rectiligne conduit à des formules analogues pour U„ et V„, et 
inversement. 

Ainsi la formule (3) 

u 2n =u n v n , 

correspond à la formule 

sin 2z = 2 sin z cos z ; 
les équations 
(6) " V n + W n = 2<f, V n -SU n = 2b% 

que l'on déduit immédiatement des formules (2) correspondent exactement 
aux relations 



cos z + V — 1 sin z = e zV ~ l , cos z — V — 1 sin z = e~ zi ~ x , 

et les formules (4) sont entièrement analogues à celles qui ont été données 
dans les Actes de Levpzick, en 1701, par Jean Bernoulli, pour le développe- 

ment de — : et de cos nz suivant les puissances du sinus et du cosinus de 

sin z 

l'arc z. Ainsi encore les formules 



(7) ' Un + &uj = 2- 1 [F m . + w nr -\ , 

que l'on déduit des relations (6) coïncident avec les formules 

(cos x + V — 1 sin x) (cos y + V — 1 sin y) =. cos (x + y) + V — 1 sin (# + #)» 

(cos a^ + V — 1 sin x) r = cos rx + V — 1 sin r#, 

qui ont été données par Moivre. 

Nous ferons encore observer que si, dans l'équation (1), on pose 
X~x\ a = a r , (3 = b r , 
les quantités a et (3 sont les racines de l'équation 
(8) X 2 =V r X— Q r . 

48 
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Par conséquent, chacune des formules qui appartiennent à la théorie présente 

peut être généralisée, en y remplaçant U n et V n par ~ et V nr , P par V r , Q par 

u r 

Q', et la différence 8 des racines a et b par la différence 8 U r des racines a et /? 
de l'équation (8). 

Les formules (4) deviennent ainsi 

Vn-lUnr _ _» F » _ 1 , « ( » — 1) (» ~ 2) A ,„ yn _ , 

?7 r ~ 1 r +_ "~ 1.2.3 "" r r 

n(n — l)(n — 2)(n — 3)(n — 4) „ 4 5 
"*" " 1.2.3.4.5 r r i "' " 



(9) 



2" - x F nr = F;' + n -^-j^ A £7? F; - 2 



w ( W _ 1) ( n — 2)(n — 3) 

+ f72":3^ A Ur K + • • • 

D'ailleurs, nous laisserons de côté, pour l'instant, les autres procédés de trans- 
formation de l'équation (1) par substitution de variable, ainsi que l'étude des 
fonctions plus générales 

AU n + BV n +C, 

dans lesquelles A, B, C désignent des nombres entiers quelconques, positifs. 
ou négatifs. 

Section III. 

Des relations de récurrence pour le calcul des valeurs des fonctions U n et V n . 

Le calcul des valeurs de U n et de V n qui correspondent aux valeurs 
entières et consécutives de n, s'effectue rapidement au moyen de formules 
entièrement analogues à celles de Thomas Simpson : 

sin (n + 2) z = 2 cos z sin [n + 1) z — sin nz , 
cos (n + 2) z — 2 cos z cos (n -\-ï) z — cos nz . 

En effet, multiplions par w n les deux membres de l'équation (1), et remplaçons 
successivement ce par a et b, nous obtenons 

a n + 2 = Pa n + i — Qa n , b n + 2 = Pb n + 1 — Qb n , 
et, par soustraction et par addition, 

U n + 2 =PU n + l -QU n , 

r n+2 = pv n+1 -Qv n . 



(10) { 
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Ces formules nous font voir que les fonctions £7" et F forment, pour les valeurs 
entières et consécutives de n, deux séries récurrentes de nombres entiers. Ces 
séries ont la même loi de formation, mais elles diffèrent par les conditions 
initiales. Nous généraliserons ces formules par l'emploi du calcul sym- 
bolique. En effet, en désignant par . F une fonction quelconque, on tire 
évidemment de l'équation (1) 

F (x 2 ) =F (Px — Q); 
si l'on remplace x par a et b, on a 

a n F (a 2 ) = a n F (Pa — Q) , b n F {V) - VF (Pb — Q) , 
et, par soustraction et par addition, on obtient les égalités symboliques 

U n F(lP) = U n F{PV— Q) , 



( U ) | y»p ( y*) = yn F (py— Q) 

dans lesquelles on remplace, après le développement, les exposants de U et de 
V par des indices, en tenant compte de l'exposant zéro. Ainsi les symboles 
U 2 et PU — Q, V 2 et PV — Q sont respectivement équivalents, et peuvent 
être remplacés l'un par l'autre dans les transformations algébriques. 

On a, par exemple, dans la série de Fibonacci, les résultats suivants 

if+* = u n -* (u + l) p , 



( '\ u n ~ p = u n (u — l) p . 

qui sont entièrement analogues à ceux que l'on peut obtenir dans la théorie 
des combinaisons ou du triangle arithmétique, et, en particulier dans la 
formule du binôme des factorielles, due à Vandermondb. 
En prenant, pour point de départ, l'équation 

on trouvera encore de nouvelles relations entre les coefficients de la même 
puissance du binôme. 

La considération de l'équation (8) conduit aux relations suivantes 

Un + ïr == 'r U n + ,. H Un j 

'n + ïr ==z 'r ' n + r H 'n ; 

qui permettent de calculer les valeurs des fonctions U n et V n qui correspondent 
à des valeurs de l'argument n en progression arithmétique de raison r. 

Inversement, on trouvera, dans la théorie des fonctions circulaires et 
hyperboliques, des formules analogues aux formules (11) et (13). 



(13) { 
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Section IV. 

Des relations des fonctions U n et V n avec les déterminants. 

On peut exprimer U„ et U„, V n et V nr au moyen de déterminants ; en 
effet, on a les formules 

f7 2 _Pf7- 1== 0, 

U,-PU 2 + QU 1 = O i 

U 4 — PU 3 + QU 2 * =0, 

U 5 -PU t + QU S * * = 0, 



U n + 1 - PU n + QU n _, 



on en déduit 



(14) U, 



n + 1 



(-1)" 



On obtient aussi 



(15) 



F.= (— 1)» 



-P, +1, 0, 0, 

+ <2, -P, +i, o, 

0, +q, — p, +1, 

0, 0, +Q, -P, 



-P, +2, 0, 0, . 

+ Q, -P, +1, 0, . 

o, +a, -p, +i, . 

0, 0, +Q, — P, . 



* =0; 



(m colonnes), 



(n colonnes), 



On vérifie les résultats que nous venons de trouver, en développant les déter- 
minants suivant les éléments de la dernière ligne ou de la dernière colonne. 

U 
Les valeurs de ~ et de V„ s'obtiennent encore au moyen des déterminants, 

en remplaçant comme à l'ordinaire P par V r , et Q par Q r . 

Enfin, nous ferons observer que ces formules sont susceptibles d'une 
grande généralisation; en effet, dans les formules (11) qui contiennent une 
fonction arbitraire, faisons n successivement égal à 1, 2, 3, .... m ; nous 
obtenons alors m équations desquelles on tirera la valeur de l'une ou l'autre 
des fonctions U et V. 
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Remarque. — On peut encore pour le développement de U n employer la 
formule suivante, 



(16) 



"n + l 



P, 


VQ, 


o, 


0, . . . . 


ve, 


P, 


VQ, 


0, . . . . 


o, 


*SQ, 


p, 


V Q, . . . . 


o, 


o, 


VG, 


P,.... 



(n colonnes), 



cependant l'emploi de la formule (14) est bien préférable. 

Section V. 

Des relations des fonctions U n et V n avec les fractions continues. 

Les fonctions U„ et V„ sont développables en fractions continues ; en effet, 
considérons l'expression 



(17) 



— — w -r t — j— 
«n h -\- a 2 






(18) { 



h 

et désignons par ^„ et S n le numérateur et le dénominateur de la n ième réduite ; 
on sait que l'on a 

-**« + 2 == Vn — 2 -tin + 1 "T ^n + 2 -"n j 

^n + 2 = "m + 2 "n + 1 "T ^ra + 2 *\ 5 

et, de plus 

(19) B n S n + 1 — B n + 1 S n = ( — l) n a l a 2 a 3 . . . . a n + 1 . 

Par conséquent, si l'on pose 

«o = h = h =....= b n = P, 

#j = $ 2 = # 3 = .... ^= $ n = — (^j 
on obtient l'expression 



(20) 



22» -M p Q_ 

TÏ7 ~ P~ 



P-Q 



P — 
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dans laquelle n désigne le nombre des, quantités égales à P. 
On a ainsi, dans la série de Fibonacci : 

( m 1 (!±_vïr + '-( î-vs )"^ 1 i__ 

2 (1 + V 5)" - (1 - V 5) n 1 + 1_ 



1 + 1 



1 + - 



dans la série de Fermât : 

(22) 



2" + 1 


— 1 




3 




2 





















2 n - 


- 1 








3 


— 2 
3 — 


2 



3 — . 



et dans la série de Pell, 

m\ (1 + V2)" + 1 -( i-V^r + 1 _ i 
k ; " (i+V2)"-(i-v"â)" " J 2-1 



2—1 

2 — 1 



2 — . 



D'ailleurs, on a généralement 

(24) -ir=" " 






1 



donc, en désignant par a la plus grande des racines, prises en valeur absolue, 
de l'équation (1), on a 

(25) Lim - n - + J = a , 

Un 

lorsque n augmente indéfiniment. Cependant, nous ferons observer que ce 
dernier résultat ne s'applique pas dans le cas des séries de troisième espèce, 
c'est-à-dire lorsque les racines de l'équation proposée (1) sont imaginaires. 

Au moyen de cette dernière formule, il est facile de calculer rapidement 
un terme de la série U„ lorsque l'on ne connait que le précédent. Soit, par 
exemple, dans la série de Fibonacci 
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u u = 7014 08733, 
= 1, 61803 39887 39894 8482 ; 



et _ 

1 + V5 



2 

si l'on calcule par les méthodes abrégées le produit «.«44, à moins d'une unité 
près, on trouve exactement, puisque u n est entier 

u a = 11349 03170. 

On peut, d'ailleurs, déterminer directement le dernier chiffre de u n ; ainsi dans 
ce cas particulier, il est facile de faire voir que deux termes, dont les rangs 
différent d'un multiple quelconque de 60, sont terminés par le même chiffre ; 
si l'on suppose alors p intérieur à 60, on peut démontrer que les derniers 
chiffres de u p et de u q sont complémentaires, lorsque la somme p + q est 
égale à 60 ; on peut donc supposer maintenant p égal à 30 ; et même p infé- 
rieur à 15, si l'on observe que les termes u 15+p et u 15 _ p ont les mêmes derniers 
chiffres, lorsque p est impair, et leurs derniers chiffres complémentaires, 
lorsque p est pair. 

On a, plus généralement, la formule 

(26) ~v^ - Vr ~ F, - or 



K — or 



v r — 



dans laquelle les V r sont en nombre n, et, lorsque n augmente indéfiniment, 

(27) Lim U(n + 1)r = or . 

A la formule (26), correspond, dans la théorie des fonctions circulaires la 
formule 



(28) 



sin (n 4- 1) z ~ 1 

1 — ' — '— = 2 cos z — - 



sin nz 2 cos z — 1 



2 cos z — 1 



2 cos z — 



dans laquelle l'expression. 2 cos z est repétée n fois. 

* Journal de Crelle, tome xvi, pag. 93 ; 1837. 
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On a aussi pour la série des V n , la relation 
(29) 







K - Q" 




n-- 






Q r 




@ 



(30) { 



dans laquelle la quantité F r est répétée n fois. 

Les nombreuses propriétés des déterminants et des fractions continues 
donnent lieu à des propriétés analogues pour les fonctions U n et V n . Ainsi 
la propriété bien connue de deux réduites consécutives, renfermée dans la 
formule (19) donne 

772 TT TJ /> - 1 

u n u n — l u n + l ^ ) 

et, plus généralement 

( fP 77 77 n(n-l)rTT2 

/qi\ 1 U «r v (n-l)ry(n + l)r H u r ) 

v* 51 / 1 v 2 V V û^ n -^ r /\TT 2 - 

\ r „ '(» — l)r r (n + l)r H ^ u rf 

on a, dans la théorie des fonctions circulaires, les formules analogues 
sin 2 * — sin (x — y) sin (x -f- y) = sin 2 y , 
cos 2 * — cos {x — y) cos (x + y) = sin 2 y , 

Il est d'ailleurs facile de vérifier immédiatement les formules (31), en 
remplaçant U, V et Q en fonction de a et b. Ainsi, on a encore 



A Vi + r = a 2 " + 2r + è 2M + 2r — 2Q n 



+ r 



A TJl = à Zn + è 2K — 2Q n ,; 
donc, par soustraction : 

a [_m +r — Q r uz] = o 2 »^ — **+•] (y — bq , 

et, par suite 

(32) U n + r — (,; U n = t/ r (-/ 2n + r ', 

on aura, par la même voie, la relation 

(33) VI +r — Q' Vl = AU r U 2n + r . 

La formule (32) donne plus particulièrement, pour r = 1, la relation 

(34) m +1 — Qm= u 2n+1 . 



